Different vacua in 2HDM 
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We discuss the extrema of the Two Higgs Doublet Model with different physical properties. We 
have found necessary and sufficient conditions for realization of the extrema with different properties 
as the vacuum state of the model. We found explicit equations for extremum energies via parameters 
of potential if it has explicitly CP conserving form. These equations allow to pick out extremum with 
lower energy - vacuum state and to look for change of extrema (phase transitions) with the variation 
of parameters of potential. Our goal is to find general picture here to apply it for description of 
early Universe. 



I. INTRODUCTION. MOTIVATION 

The Two Higgs Doublet Model (2HDM) presents 
the simplest extension of minimal scheme of Elec- 
troweak Symmetry Breaking (EWSB) allowing to in- 
clude naturally observed CP violation and Flavour 
Changing Neutral Currents [1]. Natural approach 
in its description is to derive its parameters based 
from modern data and future observations at collid- 
ers. This very approach with properties of vacuum, 
fixed by observations, is developed in many papers 
(see e.g. [2]- [4]). Such approach contains the danger that 
the discussed set of parameters allow another minimum 
of potential which is deeper than the discussed one, in 
this case obtained results correspond not real but false 
vacuum. 

The Models like 2HDM can have different vacuum 
states (minima of potential) with various physical prop- 
erties. Another approach appears - to study different 
vacua in order to confine field of possible values of pa- 
rameters of Lagrangian allowing description of data (see 
e.g. [5]-[8]). 

In this paper we consider different possible vacuum 
states in 2HDM with two goals. 

1) . We like to have a complete set of necessary and 
sufficient conditions for realization of the extrema with 
different properties as the vacuum state of the model. 
These conditions must allow to check whether the dis- 
cussed minimum of potential is global one (vacuum) or 
not. 

2) . The obtained results must allow to study (in future 
works) what can happen at variations of parameters of 
Lagrangian, related to the evolution of earlier Universe, 
as it was proposed in [9]. Let us describe this idea in 
more details. 

At first moments after Big Bang the temperature of 
the Universe T was very high, in this stage vacuum ex- 
pectation values of Higgs fields are given by minimum of 
the Gibbs potential Vq- The latter is a sum of the Higgs 
potential V(4>) and the term aT 2 4> 2 . It corresponds to the 
Higgs model with parameters varying in time. At large 
T potential has EW symmetric minimum at (</>} = 0. 
This stage describes the widely discussed phenomenon of 
inflation. 



t = o t > t, 

FIG. 1: Evolution of Gibbs potential in minimal SM 



During the inflatory expansion, the Universe becomes 
colder, and at some temperature the Gibbs potential 
transforms effectively into the well known form of the 
Higgs model with (</>) ^ - we obtain our world with 
massive particles, etc. (EWSB) (see Fig. 1). This phase 
transition determines the fate of the Universe after infla- 
tion. 

We see in 2HDM many possible vacuum states depend- 
ing on interrelation of the parameters of the potential. In 
the Gibbs potential the temperature dependent addition 
to the mass term has form T 2 aij{4>\4>j)/2 [9]. With these 
terms the interrelation mentioned above changes during 
cooling of Universe, this leads to the change of phase 
state of Universe. The sequence of phase states of Uni- 
verse during its early history, transitions among vacuum 
states with different properties can influence for current 
state of Universe [9] . 

• For these goals we consider all possible extremum 
states of the Higgs potential and after that investigate 
which of these extrema can be the vacuum state - a global 
minimum of potential. So, in this paper we discuss all 
types of extremum states in 2HDM and determine con- 
ditions when one of them is vacuum state. The obtained 
explicit form of dependence of different vacuum state en- 
ergies on parameters of Lagrangian seems to be an im- 
portant result on this way. 

In sect. II we describe the Lagrangian of the model 
and its general properties. Sect. Ill is devoted to the de- 
scription of different extrema of the potential and their 
first classification. In short sect. IV we discuss conditions 
when the electroweak symmetry point (<p±) = {<f>2} = 
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can be (local) maximum or minimum of the potential. 
In the sect. V we study the most exotic type of the ex- 
ternum - charged extremum, which is not realized in our 
world - in this extremum the interaction of gauge bosons 
with fermions will not preserve the electric charge, pho- 
ton becomes massive, etc. [7]. However it is not improb- 
able that this state was the vacuum state in some period 
after Big Bang. Then we come to the general discussion 
of "normal" neutral cxtrema in sect. VI. Their study is 
continued for the important case of the explicitly CP con- 
serving potential (with all real coefficients) - sect. VII. 
We discuss in detail uniquely defined doubly generate 
spontaneously CP violating extrema, sect. VII A and CP 
preserving extrema in sect. VII B. Then we discuss inter- 
relation among different extrema, sect. VIII. In appendix 
A we develop a special toy model, for which all the cal- 
culations can be done easily. This model gives a simple 
illustration of many general statements in the main text 
and provides an answer to problems of realizability of 
some situations. In sect. IX and X we summarize the 
results obtained and briefly discuss possible applications 
for the history of Universe. 



II. 



LAGRANGIAN 



The spontaneous electroweak symmetry breaking via 
the Higgs mechanism is described by the Lagrangian 1 



C = L s g f +C H +C Y with £ H =T-V : 



fa 



(1) 



Here L s g f describes the SU(2) x U{1) Standard Model 
interaction of gauge bosons and fermions, Cy describes 
the Yukawa interactions of fermions with Higgs scalars 
and Ch is the Higgs scalar Lagrangian; T is the Higgs 
kinetic term and V is the Higgs potential (3). In this 
paper we won't consider effects related to possible non- 
diagonal terms in T (see preliminary discussion in [4]). 

The Two Higgs Doublet Model is the simplest exten- 
sion of the minimal SM. It contains two scalar weak 
isodoublets fa- and fa with identical hypercharge. In 
particular, it is realized in MSSM. To describe Higgs po- 
tential in short form, it is useful to introduce isoscalar 
combinations of the field operators 



Xi 



0101, 



X 2 



402, 



X3 = fa\fa, X 3 * 



0201 



(2) 



The most general renormalizable Higgs potential is the 
sum of the operator — V 2 of dimension 2 and the operator 
V4 of dimension 4: 



1 Notations and main definitions follow [4], we use some equations 
from [9]. 



(3) 



V = V -V 2 (x i ) + V i (x i ); 
V 2 {x i ) = M i x i = 

- [m\ l xi+m 22 x 2 + {m\ 2 x 3 +h.c.)\ , 

Vi(xi) = AijXiXj/2 = 

X\x 2 +X 2 x 2 , . + 

= r-A 3 xia;2+A4a;3a^+ 

+ A 6 xix 3 + X 7 x 2 x 3 + h.c. 



Here i, j = 1, 2, 3, 3*, Ay = Ajj. Besides, A1-4 and 
m% are real while X5-7 and m\ 2 are generally complex. 
The field independent term Vo is added for convenience 
in future, we omit this term in many equations below. 

• The reparametrization and rephasing sym- 
metry. Our model contains two fields with identical 
quantum numbers. Therefore, the pure Higgs sector can 
be described both in terms of fields fa, {k — 1,2), used 
in Lagrangian (3), and in terms of fields fa k obtained 
from fa by a generalized rotation. The corresponding 
reparametrization symmetry was studied in [4, 6, 10]. 

In fact, in the description of reality we deal usually 
with the Yukawa sector, where right hand isosinglct 
fermion fields of each type are coupled with only one 
basic field fa or fa (Model II, like MSSM, or Model I, 
see [2]). This property becomes hidden at the general 
reparametrization transformation. The efficient form 
of the potential is that in which above property of 
Yukawa interaction is explicit. 

This efficient form of potential retains one degree of 
freedom - the independent phase transformation of fields 
fa ~ * fae iPi and corresponding phase transformation of 
parameters of potential are allowed - rephasing (RPh) 
transformation. RPh symmetry group is the sub- 
group of reparametrization symmetry group. For La- 
grangian it is one-parametric group with single parameter 
pi — p 2 . Below we have in mind mentioned efficient form 
of potential and RPh freedom for it. 

The potential with explicit CP conservation is that 
with all real Xi, m\ 2 (or with parameters which can be 
transformed to real ones by single RPh transformation). 

• The results for the most general Lagrangian (pre- 
sented below) often have very complex form. Main fea- 
tures of physical picture are seen in more simple potential 
with softly broken Z 2 symmetry 2 , in which Xq = A 7 = 0. 
We will discuss many results using final equations for this 
very potential, often - in the explicit CP conserving case. 

• Positivity constraints. To have a stable vacuum, 
the potential must be positive at large quasi-classical val- 



2 For the most general Higgs potential loop corrections with Ag, A7 
generally mix <pi and <f>2 fields even at small distances, it results 
in breaking of the mentioned Model II or Model I form of Yukawa 
interaction. This breaking is absent for the potential with softly 
broken Z2 symmetry, in which the kinetic term has diagonal form. 
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ues of fields \<pk\ (positivity constraints) for an arbitrary 
direction in the (fa, fa) plane. It means that upon re- 
placement of the operators Xi with some numbers Zi 

Vi(zi) > if z\ , z 2 > , ziz 2 - z 3 z 3 > . (4a) 

This condition limits possible values of A^. For the po- 
tential with softly broken Z 2 symmetry (A6 = A7 = 0) 
such limitations have form (see e.g. [11, 12]) 



Ai>0, A 2 > 0, \A1A2 + A 3 > 0, 
v / X^ + A 3 + A 4 - |A 5 | >0. 



(4b) 



III. EXTREMA OF POTENTIAL 

The extrema of the potential define the values (fa ,2} 
of the fields <f>i^ 2 via equations: 

dv/dfa\^ =i4>i) = o, dv/dcfil^^) = o. (5) 

These equations have the electroweak symmetry conserv- 
ing (EWc) solution (fa) = and the electroweak symme- 
try breaking (EWSB) solutions. Here and below nota- 
tion (F) mean numerical value of the operator F at ex- 
ternum. In general, there are many EWSB extrema. We 
label these extrema by an additional subscript, if neces- 
sary; e.g. (F) n means value of (F) in iV-th extremum. 

We consider also the values yi of operators Xi at the 
extremum points. In the discussed tree approximation 
(mean field in the statistical physics) we have 



Vi,N 



(Xi) N = (<j> a % (fa) 



)N 



for 



(7) 



In each extremum point these values obey inequalities 
following from definition and Cauchy inequality, written 
for important auxiliary quantity Z: 

yi > , y 2 > , Z = y x y 2 - y 3 y 3 > . (6) 

• Classification of EWSB extrema. It is useful to 
define quantities 

T a = (dV/dx a ) = A aiVi - M a (a = 1, 2, 3, 3*) , 

Ti, 2 are real, T 3 , = T* . 

In these terms eq. (5) have form 

(dV/dfa\) = (fa)T 1 + (fa)T 3 = 0, 

(dV/d4)= (fa)T 2 + (fa)T 3 , =0, 

(dV/dfa) = (fa^T, + (fa)^T 3 , =0, 

(dV/dfa) = (fa^T 2 + (fa)^T 3 = 0. 

These equations can be easily transformed to equations 
for yi. For example, 

(fa)^(dV/dfa\) = yiT 1 +y 3 T 3 = 0, 

(fa)^(dV/dfa\)=y* 3 T 1 +y 2 T 3 = Q, 

(fa)t(dV/dfa\) = y 2 T 2 + y* 3 T 3 , = , 

(fa)t(dV/dfa\) = y 3 T 2 + Vl T 3 , = . 



(8) 



One can consider each pair of these equations as a sys- 
tem for calculation of quantities Ti via yi. The deter- 
minant of these systems are precisely Z = y\y 2 — y 3 y 3 - 
Therefore, it is natural to distinguish two types of ex- 
trema, with ("charged extrema" with Ti = 0) and 
with Z = ("neutral extrema" with T 4 ^ 0). 

• For each EWSB extremum one can choose the z 

axis in the weak isospin space so that (fa) = ^ ^ with 

real v\ > (choose "neutral direction"). In such basis 
(fa) has generally an arbitrary form. Then, after this 
choice the most general electroweak symmetry violating 
solution of (5) can be written in a form with real v\ and 
complex v 2 : 



^(O' ^ V2W (9) 
with vi — \vi\, v 2 — 1^2 , 



without loss of generality we can consider only real pos- 
itive u. 

At u 7^ we have Z > - charged extremum, at u = 
we have Z = - neutral extremum. 

• Mass matrix at each extremum is given by the 
decomposition of the fields fa near the extremum point. 
Since the position of an extremum point relative to origin 
of coordinates selects some direction in the isospace, this 
matrix has different elements for different components of 
fields 



Mi 



d 2 v 



ij,af3 



b?dfa] a 



,(i,j =1, 2, a,p = +, 0). (10a) 



Direct differentiation gives 

M lhaf3 = T 1 6 a p + Au(fa0)Hfa) a + 
+A* 13 (fae)Hfa a ) +A 3 i(fa f 3)Hfa a )+ 
+A?,y(fap)Hfa a ) , 

M l2 . afi = T 3 S aP + Ai 2 (fa^(fa a ) + 

> + 



»2a 



(10b) 



+Al3*(fap)Hfa a ) +A 32 (02/3) 
+A33* (fafj)Hfa a ) , 

M 2 l, a p = T 3 ,5 af ) + Ai 2 (02/3)(^la> t + 
+Ai3*(01 /3 )W>la) t +A 3 . 2 (02/3)(02 Q ) t + 
+A33- (fap) (fa a ) f , 

M 22Mi = T 2 S a p + A 22 (fa ) f (faa ) + 

+A 23 (fap) J *(fa a ) +A 3 * 2 (fa^(fa a )+ 

+A33- (fa f i)Hfa a ). 

• The distances from some extremum and be- 
tween two extrema are defined as 

V(faN) = 

= {fa(fa)N-fa(fa)N)' { (fa (fa) N - fa (fa) n)) = 

= x x y 2 + x 2 yi - x 3 y 3 * - x\y 3 , (11) 
V(I, II) = (V(fa II)) j = (V(fa I)) n = V(II, I) = 
= \((faj)(foji)~(faj)(faji))\ 2 - 



4 



Note that T>(fa N) = at fa = 
for a pair of different extrema. 
• The extremum energy is 

= vm N ) = 

= V(y itN ) = V - V 2 {yi, N ) + V 4 ( yi , N ) . 



(fa)N and T>(I,II) > quantities yi in a following way 



(12) 



According to theorem on homogeneous functions, in 
each extremum point 

^ijyi,NVj,N = Miy itN , or equivalently 

V 2 ((fa) N ) = m((<f>i) N ) => (13) 
=> £ e N xt =V - V 4 (y hN ) = V - V 2 ( yi ,N)/2 . 

The extremum with the lowest value of energy (the 
global minimum of potential) realizes the vacuum state 
of the model. Other extrema can be either saddle points 
or maxima or local minima of the potential. It can be 
established, in particular, by the study of the effective 
mass matrix at these extrema. 

• Decomposition around EWSB extremum. 
Our potential can be rewritten as a sum of extremum 
energy and terms vanishing in the extremum point to- 
gether with their derivatives in fa . This polynomial with 
terms up to the second order in X, can be written as a 
sum of polynomials of second and first orders in Xj van- 
ishing together with their derivatives in the extremum 
point. The form of second order polynomial is fixed by 
a quartic terms of potential, it can be only V±{xi — y{). 
The residuary first order polynomial in Xj must be pro- 
portional to T>(fa N). Therefore 



V 



£ e N xt + V 4 (xi - Vi, N ) + K ■ V(fa N) . (14a) 



It means that the mass terms of (3) can be written via 



J 



m ii = 2 ( A U/i + hVz + A 62/3 + Agi/3* - yiK) , 
m\ 2 = 2 (A 3 yi + \ 2 y 2 + A 7 y 3 + A^ 3 « - y{R) , (14b) 



rn 



12 



2 (A 6 j/i + A 7 y 2 + A 5 7/ 3 + A 4 y 3 « - yyK) . 



The differentiation of (14a) gives for Tf. 
Ti = y 2 TZ , T 2 = yi TZ , T 3 - -y* 3 TZ , T 3 . - -y 3 TZ . 

In accordance with (8) for charged extremum we have 
from here 1Z — 0. 

For the neutral extremum the Higgs fields mass matrix 
(10b) for the upper (±) components can be written as 



M 



++ 



Ti T 3 

T 3 , T 2 



- / 2/2 



"2/3 



-j/3 yi 



K. 



At Z — determinant of this matrix equals to 0. There- 
fore, one eigenstate of this matrix equals to 0. It de- 
scribes massless combination of charged Higgs fields (well 
known Goldstone state). The second eigenstate of above 
matrix describes the physical charged Higgs boson with 
mass 



Tr M++ =T 1 +T 2 = ( Vl + y 2 )K . 



This quantity is positive for the minimum of the poten- 
tial, it can be negative in other extremes. Finally, we 
obtain 



K 



n = o 

M% ± 



2/i + 2/2 



for charged extremum 
for neutral extremum N . 



(14c) 



N 



IV. EW SYMMETRY CONSERVING (EWC) POINT 



The EWc point (fa) = (fa) 



is extremum of potential. Depending on mf • it has different nature: 



At 
At 



m\ x < , m\ 2 < , 
m[ x > , m 22 > 



and rn\ x rw 22 > \m\ 2 



2 12 



According to [10] no other extremum can be a maximum of potential. 



— minimum, 

— maximum , 

saddle point. 



(15) 



V. CHARGED EXTREMUM 

We consider now the extremum which appears at 

Z = 2/12/2 - 2/3*2/3 ^ => u^Q. (16) 

If this extremum realizes the vacuum, it is not possi- 
ble to split the gauge boson mass matrix into the neu- 
tral and charged sectors, the interaction of gauge bosons 



with fermions will not preserve electric charge, photon be- 
comes massive, etc. [7]. That is the reason why this ex- 
tremum is called the charged extremum. We label quan- 
tities related to this extremum by a subscript N = ch, if 
necessary. 

Certainly, this case is not realized in our World. Nev- 
ertheless, it is interesting to consider main features of the 
case when this extremum is the vacuum state in respect 
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to the opportunity of different scenarios in the Early Uni- 
verse. 

In the considered case eqs. (8) for the extremum of the 
potential have form 



T x = X lVl + X 3 y 2 + X* 6 y* 3 + X 6 y 3 - ro^/2 = 0, 

T 2 = X 2 y 2 + X 3 yi + X 7 y 3 + X 7 y 3 - m\ 2 j2 = 0, 

T 3 = A 4 2/3+A 5 y 3 + A62/i + A 7 y 2 - m\ 2 /2 = 0, 
T y = A 4? y 3 + A^3* + A^i + X* 7 y 2 - mf 2 /2 = 0. 



(17) 



That is a system of linear equations for y^. It can have 
only one unique solution (except some degenerate cases) 

system can have only one charged extremum. 

This system has solution at arbitrary parameters of the 
potential. However, in accordance with (9), it describes 
an extremum of the original potential (3) (we define some 
quantities Vi and u) if only the obtained values yi i2 obey 
inequalities y\ > 0, y 2 > 0, Z > (6). In this case 



u=^2Z/y~ u e^=y 3 /y* 3 . 



(18) 



Inequalities (6) determine the range of possible values of 
Aj and m?- where the charged extremum can exist. 

• Condition for minimum. In the discussed case 
the potential (3) can be rewritten in the form (14) with 
1Z = 0. The charged extremum is minimum of the po- 
tential if the quadratic form V±{xi — yi^h) is positively 
defined at each classical value of operators Xi. This con- 
dition differs from the positivity constraint (4a), since 
quantities Zi — Xi — yi^u do not need to satisfy condi- 
tions given in this constraint. Here z\ and z 2 are some 
real quantities (positive or negative) and z 3 — z 3 , is an 
independent complex quantity. Therefore, the condition 
for the charged minimum is (see e.g. [10]) 

Vi(zi) > at arbitrary real Z\, z 2 , complex z 3 . (19) 

• The case of softly broken symmetry 

(A 6 = A7 = 0). The main features of the solution 
are seen in the case of the soft Z 2 symmetry violation. 
In this case solution of equation (17) has form 



yi 



m\ 1 \ 2 — m\ 2 \ 3 — m\ x X 3 +m\ 2 X\ 

2(AiA 2 -A2) ' V2= 2(X 1 X 2 -X 2 3 ) 



1/3 = 



ml 2 X* 5 +m*^X 



(20) 



12 A 4 



2(A| - A5A5) 



Conditions (6) limit the domain in the space of param- 
eters, where the charged extremum can exist. 
The first condition in (6) reads as 



(mf 1 A 2 - m^ 2 A 3 )(m| 2 Ai - mf 1 A 3 ) 



(A1A2 



A|) 2 



m 



12 A4 



12 ^5 



A4 — A5A 5 



(21a) 



> 0. 



The specific form of condition y\ 2 > depends on 
the value of parameter A 3 admissible by the positivity 
constraint (4b): 



m\ x X 2 >m\ 2 X 3 , 
m 22 Ai >to 2 1 A 3 
m\ x X 2 <m\ 2 X 3 , 



at 



at 



I A3 



\A1A2 
A 3 



m 22 Ai<mf 1 A3 \/AiA 2 



< 1: 



> 1. 



(21b) 



If A5 and m\ 2 are real (explicitly CP conserved poten- 
tial), condition (21a) forbids also small values of A 4 + A 5 . 

The extremum energy in this case is subdivided into 
the sum of Z 2 symmetry conserving term and Z 2 sym- 
metry violating term: 



cext 

C ch — 



m\ x X 2 +m\ 2 X\ — 2m\ 1 m\ 2 X 3 
8(A!A 2 -A2) 
2m^ 2 m\ 2 X4 L — m\ 2 X5 — m \ 2 ^\ 
8(A| - A 5 A*) ' 



(22) 



This extremum is the minimum of the potential 
(charged minimum) if condition (19) is satisfied. In the 
considered case it is easy to obtain (with argumentation 
similar to that in [12]) that this condition is satisfied if 

Ai>0, A 2 >0, v / AiA 2 + A 3 > 0, A 4 >|A 5 |. 

Note that latter condition guarantees negativity of sec- 
ond item in (22). 



VI. NEUTRAL EXTREMA, GENERAL CASE 

Other solutions of the extremum condition (5) obey a 
condition for U(l) symmetry of electromagnetism, that 
is solution with 



Z = ym - y 3 y 3 = 



2}: 







V2 \ v iJ ' V2 \ v 2 = M e *V ' 

another parameterization: (23) 
v\ — v cos(3 , v 2 — v sin[3 , 
^ = 2(y 1+ y 2 ) = 2(|<«M| 2 + K<M| 2 ) >0, 
e 2 * = 2/3/2/3* • 

In this case quantities y, are not independent. There- 
fore, the field values at the extremum point cannot be 
obtained by minimization of form (12) in j/j. In these 
terms system of equations for v.e.v.'s has form (8) with 
solutions yi = yi_ n - It is important to note that the num- 
ber of independent parameters here is 3 (not 4). Those 
are the real quantities y\, y 2 and the phase difference 
of values of fields at the extremum point £ (not 
separate phases of these values!). 
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• Charged Higgs mass. It is useful to reproduce 
here the equation for the charged Higgs mass, given e.g. 
in [4] (cq. (4.3)). 



Re([m 2 12 -2(X 6 y 1 + X 7 y 2 )] e 
A 4 + Re(X 5 e- 2 ^) 



■<0 



(24a) 



In addition, we present corrected eq. (4.5c) from [4] 
which gives mass of either pseudoscalar Higgs in the case 
of CP conservation or intermediate quantity obtained at 
partial diagonalization of mass matrix in the case of CP 
violation 



M 3 



Mi 



M 2 H± 



v 

T 



(A 4 -i?e(A 5 e- 2l «)) . (24b) 



• For the Higgs potential of general form we have no 
idea about classification of neutral extrema. However, 
if CP conserving extremum (with no scalar-pseudoscalar 
mixing) exists, there is a basis in (<pi, r/> 2 ) space in which 
potential has explicitly CP conserving form [5] , [4] (with 
all real A i; to?-). Using such a form of potential, the 
subsequent useful classification can be introduced, in this 
important case. 



VII. NEUTRAL EXTREMA, CASE OF 
EXPLICIT CP CONSERVATION (REAL A;, m? ) 



parameters of the potential become complex, giving CP 
violation in the Higgs sector (mixing of scalar and pseu- 
doscalar neutral Higgs bosons) - sec for details e.g. [4]. 
That is the reason why this extremum is called the spon- 
taneously CP violating (sCPv) extremum [1, 7]. 
The substitution of cos£ (26) into (25) results in 



E. 



sCPv 



My 2 + A 2 y| 



m^yi + m\ 2 y 2 



F A3452/12/2 - 
2(A 6 yi + X 7 y 2 )} 2 



(27) 



(here A 345 



8A 5 

A3 + A 4 — A5) . 



Further minimization gives the system of linear equa- 
tions dE s cpv/dyi — 0, dE s cp v /dy 2 = (we don't 
present this general system only due to its bulkincss). 
Only one solution of this system exists, so yi, and 
cos £ are described by parameters of the potential 
unambiguously. (This conclusion can be obtained also 
from the description of [10].) 

The discussed extremum can be realized only in the 
range of parameters of the potential obeying inequalities 



'12 



2(A 6 j/i + X 7 y 2 ) 



4A5V2/12/2 

yi > o. 



= icosei < 1, 



(28) 



y 2 > 0. 



The change (cf> 2 ) -> (fa)* (£ - 
the extremum energy (25). Therefore 



£) does not modify 



In accordance with definitions (9), we have for each 
solution J/3 = ^yiy 2 1 - ■ In the discussed case the ex- 
tremum energy (12) is transformed to the form 

1 



£ ext = -- {ml 1 y 1 + ml 2 y 2 + 2ml 2 ^yry 2 ~cos£} + 



Ai 



+yy? + yS/i + ( A 3 + A 4 )?yiy 2 + 
+A 5 yiy2Cos2^ + 2 (A 6 j/i + A 7 y 2 ) ^/yiy 2 cos£,. 



(25) 



Now we find extrema in coordinates y\, y 2 , £. We start 
from the minimization in £ at fixed yi. It gives two types 
of solutions: 



[A] : cos£ 



rn 



12 



2(A 6 j/i + X 7 y 2 ) 



4A 5% /yIy 2 " 
[B] : sin£ = 0. 



(26) 



For the discussed explicitly CP conserving potential this 
equation is equivalent to the constraint eq. (3.11) ob- 
tained in [4]. 



A. Spontaneously CP violating extremum 



if <j)\ — (4>\), 4> 2 = (4> 2 ) is the extremum 
of potential — (<j>\), 4> 2 = (^2)* is also 

the extremum; these two extrema are de- 
generate in energy [1] and define two "direc- 
tions" of CP violation, "left" and "right", with 



(29) 



c , (m\ 2 - 2(A 6 yi + X 7 y 2 )\ 

t = ± arcc °s tt — ■== - 30 

V ^^/yrm J 

• Our potential (25) is a second order polynomial in 
cos £. The sCPv extremum (if it exist) can be a minimum 
only if A5 > 0, in accordance with [6]. 

• Substituting the cos£ (26) into (24a) we obtain al- 
ternative form for the mass of and M 33 in sCPv 
extremum 3 



M H±,sCPv = y( A 5 - A 4) : 

-M33, sCPv = f 2 A 5 sin 2 £ . 



(31) 



To realize minimum of potential, all mass squared eigen- 
values and all diagonal terms of squared mass matrix 



The extremum point (26 [A]) describes a solution with 
complex value of field v 2 at real parameters of the poten- 
tial. The rephasing transformation cf> 2 — > 4> 2 e~ 1 ^ trans- 
forms the potential to the real vacuum form, in which 



3 The same equation for the mass of was found in [7] and (for 

the case of exact Z2 symmetry) in [14]. Note that in the discussion 

of sCPv extremum authors of [15] used equation for M\, which is 
incorrect in this case. 
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must be positive. Therefore, necessary conditions for re- 
alization of sCPv minimum arc 



A 5 > , A 5 > A 4 . 



(32) 



Note that the mass matrix for neutral Higgses in this 
scPv extremum can be written as (sign • • • mean the 
diagonal symmetric quantity from the upper right corner 
of matrix) 



extremum (28) can be written in the form 

ml 1 \2 > A 3 4 5 ra 22 
m\ 2 X\ > A 345 to 2 i 




at \7AiA 2 > |A 345 | , 



mf 1 A 2 < A 3 4 5 to| 2 
m\ 2 X\ < A 34 5m? 1 



(35) 




w 2 (X 5 si3 cos £ + AeC/3) sin £\ 
v 2 {X 5 Cfj cos£ + A7S/3) sin£ 
A 5 sin 2 £ 



— 2 < \ t~ a t V Ai A 2 < A 3 45 , 



in 



22 



m\ 2 < 16A^i2/2 • 



It shows that - as it is naturally expected - the CP- 
violating field mixing is weak if the phase £ is small. 
Moreover, in the case of soft Z2 symmetry violation 
(A 6 = A 7 = 0) physical CP violation is absent at £ = n/2 
(that corresponds to the exact Z2 symmetric case, with 
m\ 2 = 0, A 6 = A 7 = 0) 4 . 

a. The case of softly broken Z 2 symmetry (\q — 0, 
A 7 = 0). The analysis becomes more transparent for 
the softly Z 2 violating case at real A 5 , m\ 2 . The cor- 
responding linear system is \\y\ + X 345 y 2 = 
A22/2 + A 3 45j/i = m 22 /2 with solution 



B. CP conserving extrema 

The solution (26 [B]) describes extrema that corre- 
spond to £ = and £ = tt. The case £ = 7r can be 
obtained from the case £ = if we allow V2 to be nega- 
tive, i.e. allow tan to be negative. Therefore, without 
loss of generality we consider below the only case with 

e = o. 

In these cases CP violation does not appear (CP con- 
serving — CPc — extrema). The extremum condition 
(8), written for Vi = \f2yi, has form of the system of two 
cubic equations: 



ml 1 X 2 ~X 3i5 ml 2 
2A. S 

=> tan 2 = 



V2- 



2A S 



■2 , _ TO22A1 - \zA5 m \\ 

mf 1 A 2 - A 34 5to| 2 



(33) 



where A 
cos£ 



A1A2 — A 



345 > 



'12 



4A 5 Vyi2/2 



In this case the extremum energy is obtained easily from 
(12) (cf. cq. (22)): 



m'l 1 \2+m^ 2 X 1 ~2ml 1 ml 2 X 3 45 mf 2 
' sCPv = 8A S 8AT (M) 



mlivi 



m\ 2 V2 



- m\ 2 v 2 = Xivf + X345V1VI+ 
+ (3X 6 v 2 + X 7 v 2 )v 2 , 



X2V 2 + A 3 45^ 2 W2 + 

3A 7 w|)wi , 



rn\ 2 v\ ■ 
+(A 6 ^ 2 



A 3 45 — A 3 + A4 + A5 



(36) 



Rewriting this system with parametrization 
vi — vcos0, V2 — wsin/3, we express the quantity 
v 2 via t = tan/3 and obtain the equation for t, similar to 
equations presented in [8] : 



:(1+* 2 ) 

(1+t 2 ) 



tmi 



Ai + A 345 i 2 + 3A 6 t + X 7 t 3 



trrinn 



m\ 2 



X 2 t 3 + A 345 i + A 6 + 3A 7 t 2 



(37a) 



Taking into account the positivity constraint (4b), the 
necessary conditions for realization of the CP violating 



4 In the exact Z2 symmetric case sign of A5 can be changed by sim- 
ple rephasing without change of other parameters of potential. To- 
gether with condition (32) it means that formal sCPv extremum in 
this case give no CP violation. 



and 



(A 2 m? 2 -A 7 m 2 ! 2 )i 4 + 
+(A 2 m 2 1 + 2A 7 mf 2 - X 3Ab m 22 )t 3 + 

+3(A 7 m 2 1 -A 6 TO 2 . 2 )t 2 + 
+ (A 345 m 2 1 - 2A 6 TO 2 2 - X 1 ml 2 )t+ 

+ (A 6 mf 1 - Xim\ 2 ) = . 



(37b) 
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It is easy to obtain that in this parametrization 

.,2 v 2 t 2 



v{ = 2 Vl = 



4 



1 + t 2 ' 
V1V2 = 2y 3 = 2yl = 



2j/2 = 

v 2 t 

l + t 2 



l + t 2 '' 



(37c) 



Now one can rewrite extremum energy (13) for discussed 
cxtrcmum in the form 



(raj! + tm\ 2 ) {fn\i + 2tm\ 2 + t 2 i 
£cPc = 8(A X + A 345 t 2 + 3A 6 t + A 7 t3) 



22 , 



(38) 



By construction, one should consider only real solu- 
tions of equation (37b) satisfying v 2 > 0. The equation 
of fourth degree might have 0, 2 or 4 real solutions (in- 
cluding accidental degeneracy) . Taking into account pos- 
sible negative values of v 2 given by (37a) one can state 
carefully that there could be up to 4 CPc extrema. 
If necessary, we label different extrema of such type with 
an additional subscript N = I, II, .... 



• Note that at 



ii 



'22- 



Ai 



•^2, A 6 = A7 our po- 
tential has additional symmetry fa <-> fa. The potential 
has extrema keeping this very symmetry (t = ±1) and 
those where this symmetry is spontaneously broken. The 
latter states are degenerated in energy like sCPv states 
discussed in sect. VII A. If these states form minimum of 
potential, we deal with two degenerated minima (cf. [10]). 

Weak violation of mentioned symmetry destroys dis- 
cussed degeneracy. Therefore, in the general case in some 
region of parameters our system can have two CPc min- 
ima simultaneously. 

• In the case of soft Z2 symmetry violation 
(A6 = A7 = 0) equations (37) become: 



= {l+t 2 ) 



tra 



A345^ 2 + Ai 



A 2 m 2 2 i 4 + (Asm 2 ! - A 345 m| 2 )i 3 + 
+ (A 345 m 2 1 - Ximl 2 )t - Aito 2 2 = . 



(39a) 



(39b) 



VIII. VACUUM AND OTHER EXTREMA 

In this section we calculate the difference of extremum 
energies for different extrema. For this purpose we ex- 
press extremum energy in the extremum II via parame- 
ters of extremum I and vice versa with (14): 

£ ext = £ ext + y 4 ( yul _ ^ + Kj . V( j^ J) _ 

EY* = £ff + V 4 (y u - y UI ) + K H ■ V(I, II) . 

The quadratic polynomial Vi{zi) = Vi{— z{) (simultane- 
ous change of signs of its arguments). The distance be- 
tween two extrema T>(I, II) is symmetric and positive by 



definition (11). Therefore by substraction of one equa- 
tion from another we obtain 5 



A£(II,I) 



cext 
C II 



ext 



Hn)-V(I,II) (40) 



with 1Z = for charged extremum and TZ = M H ±/v 2 \ N 
for neutral extrema (14). Therefore, the quantity 
(1Zi — TZu) determines the hierarchy of extrema. 
• EWc = (02 ) = 0) and EWSB extrema. 

Rewriting eq. (14) for any EWSB extremum one can 
obtain 

= SeWc = Sew SB + Vi{yi,EWSB ~ y%,EWc) 



with yi,EWc — => 
£ewsb 



EWSB 



)<0. 



Therefore, 

1 . If EWSB extremum is a minimum of potential, then 
the EWc extremum has higher energy. 

2. If the EWc extremum ((<j>i) = (fa) — 0) realizes the 
vacuum state (it can happen only at m 2 1; m\ 2 < 0) 
all EWSB extrema are either saddle points or local 
maxima of potential. 

• Neutral extremum and charged extremum. 

1. According to eq. (14), the Higgs potential can be 
written as a sum of charged extremum energy £*j^ 
and operator V^(xi — yi.ch)- The latter is a poly- 
nomial of the second degree in Zi — Xi — yi :C h- If 
charged extremum is a minimum of the potential, 
the quantity V±(zi) is positively definite at arbi- 
trary real z\ y2 and complex Z3 (19). Since V4 is 
quadratic form in z%, this quantity is positive also 
in the points, correspondent all other extrema of 
potential. Therefore, if charged extremum is 
a minimum of the potential, it is the global 
minimum the vacuum state 6 . 

2. Besides, the difference of the cxtrcmum energies for 
neutral and charged extrema (40) can be written 
as A£(ch,n) = (M 2 H± n jv 2 n )V{ch,n) (see also [7], 
[10]). Therefore: 

2a. If the neutral extremum is a minimum of po- 
tential, i.e. Mfj ± > 0, the charged extremum has 
a higher energy [7]. 

2b. If the charged extremum realizes the mini- 
mum of the potential, all neutral extrema are sad- 
dle points or local maximums (since it can take 
place only at M\ ± < for all n). 

• Two neutral extrema I and II. 

For two neutral extrema the difference of extremum 
energies (40) is 



A£(II,I) = ( 



■V(II,I). (41) 



5 In particular cases the same kind of equations was obtained in [7]. 

6 This conclusion can be obtained also from discussion of ref. [13]. 
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Therefore, in particular, 

1. If extremum I is minimum (Mfj ± T > 0) and ex- 
tremum II is not a minimum with ^ < 0, 
extremum II is higher than extremum I. 

If two minima of the potential exist with energies 
Ex and £zz, the energy interval (£x, Six) cannot 
contain saddle points or maxima. 

2. For two neutral minima of potential or a minimum 
and a saddle point with Mfj ± N > 0, the deeper 
(a candidate for the global minimum - the vac- 
uum) is the extremum with the larger value of ratio 

M h ± , n / v %- 

• Explicitly CP conserving potential. Neutral 
extrema. 

1. sCPv and CPc extrema. 

We have M 2 H± /v 2 \ sCPv = (A 5 - A 4 )/2 (31) . In 
accordance with (24b) we have for CPc extremum 
M 2 H± /v 2 \ CPc -{\^-\ i )/2 = M 2 A /v 2 \ CPc . It allows 
us to rewrite eq. (41) in the form (see [7]) 

A£ sCPv ,cpc = Ml/v 2 \ CPc ■ V(sCPv, CPc) (42) 
with positive V(sCPv, CPc). 

Therefore, similarly to the comparison of charged 
and neutral extrema: 

la) If system has a sCPv minimum, i. e. 
M^/v 2 \cpc < 0, the minimum is the vacuum. All 
CPc extrema are saddle points, not minima, 
lb) If system has a CPc minimum i. e. 
M\/v 2 \cp c > 0, the sCPv extremum cannot be 
a minimum, it is a saddle point. 
The toy model gives good illustration for this state- 
ment (see Fig. 2). 

2. Two CPc extrema I and II in the case of the 
softly Z 2 violating potential. 

In this case one can use the eq. (24a) 
2Mf {± /v 2 — (m\ 2 / 'viv 2 ) — A4 — A 5 . It allows 
to transform the mass factor Aiuj in (41) to the 
form 



M 



11, 1 



2l>l,/t>2,/ 2^1/7^2,77 



(43) 



or - with the aid of (37) - to express Ainj via 
solutions of equation for tan (3 = t. 
So that, in this case the hierarchy of extrema can 
be established without direct calculation of M^ ± . 



IX. GENERAL PICTURE 

Let us summarize the picture obtained. 
In our discussion of a separate extremum we have in 
mind such a choice of the z axis in the weak isospin space 

that in this extremum (fa) = ^ with real v\ > (9). 

We think that it is unnecessary to search for the basic 
independent form of the results. In our opinion there is 



a efficient form of the potential, determined by the form 
of the Yukawa sector (see discussion after (3)). 

There are two very different types of extremum of the 
potential in 2HDM (9) - the charged extremum with 
u 7^ (it docs not describe the modern reality) and the 
neutral extremum with u = 0. 

A. General case. 

1) The Electroweak symmetry conserving extremum 

(((fi) — ((P2) = 0) can realizes the vacuum state if 

TO ii> m 22 < u ; I m i2l 2 < m ii m 22- I n tn i s case au other 
extrema are saddle points. 

2) The charged extremum is determined by parame- 
ters of the potential uniquely by the system of linear 
equations (17), (18). It exists if solutions of this system 
obey conditions (16). 

If the charged extremum realizes the minimum of the 
potential, it describes the global minimum - the vacuum 
state. In this case all neutral EWSB extrema are saddle 
points of the potential. 

3) The number of neutral extrema is more than one 
(in addition to EW symmetry conserving extremum 

(fa) = (fa) = o). 

The potential can have simultaneously two neutral ex- 
trema I and II. If M 2 H± Jvj - M 2 H± n /v 2 II > 0, the state 

1 is below state II (state II cannot be vacuum) (41). 

B. Explicitly CP conserving potential (with all 
real coefficients in (3)). In this case one can distinguish 
a CP conserving (CPc) extremum with zero phase differ- 
ence between the values (fa) at the extremum point and 
spontaneously CP violating (sCPv) extrema, in which 
the phase difference between the values (fa) is nonzero, 
the latter generates neutral Higgs states without definite 
CP parity. Total number of extrema in this case can be 
up to 8 (0 or 1 charged extremum, up to 4 CPc extrema, 

2 or sCPv extrema, 1 EWc extremum). 

1) The sCPv extremum is determined by the parame- 
ters of the potential uniquely by system of linear equa- 
tions. It exists if solutions of this system obey the condi- 
tion (28). The sCPv extremum state is doubly degener- 
ate in sign of phase difference between the values of fields 
at the extremum point. 

For this extremum to be minimum it is necessary to 
have A5 > and A5 > A4. If sCPv extremum realize 
minimum of potential, it is the vacuum state (doubly de- 
generated). In this case all EWSB extrema are saddle 
points. (In MSSM with loop correction A 5 < 0, there- 
fore the sCPv extremum cannot be the vacuum (see e.g. 
[15]).) 

2) System can have more than one CPc local minima, 
the vacuum state is the lowest among them. For the 
important case of softly broken Z 2 symmetry, if in this 
case we have two CPc minima of the potential, I and 
II, the eq. (43) means that the state I is below state II 
and can describe vacuum if 

m\ 2 l(v\sin2^j) - m\ 2 j (v 2 n sin2(3 u ) > . (44) 

C. The decomposition of potential near extrema (14) 
seems to be useful for phenomenological analysis. 
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D. For explicitly CP conserving potential we have 
found explicit equations for cxtrcmum energies (22), (34), 
(38) via parameters of potential and set of necessary con- 
ditions for realization of charged or sCPv extrema. These 
equations allow to pick out extremum with lower energy 
- vacuum state and to look for phase transitions at the 
variation of parameters of potential. 

X. WHAT NEXT 

Now we discuss briefly possible effects of radiative cor- 
rections and present first ideas about and the impact of 
our analysis on cosmology. These are the problems for 
future studies. 

• Radiative corrections. With radiative (loop) cor- 
rections main qualitative features of obtained picture 
will be changed weakly (provided these corrections arc 
small 7 ). These corrections are important if they vio- 
late some artificial symmetry of the potential. In our 
case that is explicitly CP conserving form of poten- 
tial. Radiative corrections contain contributions e.g. 
of light quarks, having imaginary parts for the consid- 
ered mass interval. (The simplest example gives cor- 
rection to A 5 term, obliged by interaction with 6-quark. 
Very rough estimate of loop correction gives additional 
ImXs < (mb/v) A {mb/Mh) 2 ~ 10~ 10 , where factors m^jv 
are from Yukawa coupling and factor (rrib/Mh) 2 - from 
loop integral itself.) These imaginary parts eliminate de- 
generacy of the sCPv extrema in accordance with the 
arrow of time, and it is natural to expect that the en- 
ergy difference between these two states is small - we 
deal with almost degenerate states. In simple words, one 
can write that the phase with left violation of CP is real 
vacuum. 

• Possible effects for cosmology. 

1) Temperature dependence of the Gibbs po- 
tential is determined by standard methods of statisti- 
cal physics (see e.g. [16]). In the first approximation of 
perturbation theory only mass term modifies. At large 
enough temperature T we have Am?- = a^T 2 (<j)\fa) /2. 
Higher order corrections modify these and change 
weakly (no more than logaritmically) parameters \. 
Therefore general features of phase transitions during 
evolution of Universe can be analyzed in terms of vari- 
ation of mass term of potential (3) with the aid of eq-s 
(22), (34), (38) (sec [17]). This picture can be modified 
by variation of form of Gibbs potential near phase tran- 
sition. 

Moreover, in vicinity of phase transition all processes 
become slower, and adiabatic approximation for calcula- 

7 In the discussion of Higgs sector in MSSM 1-loop corrections are 
not small. 



tion of the Gibbs potential become invalid (evolution of 
Universe can be faster than transition to the thcrmody- 
namical equilibrium). 

2) If during evolution, Universe changes different 
phases of 2HDM, these phase transitions look as tran- 
sitions of the second order as long as we consider effects 
in the tree approximation (fluctuations = multiloop ef- 
fects can change type of transition). In each phase during 
cooling the masses of the particles of matter (which are 
given by values (fa) and (fa)) evolve and even their hi- 
erarchy can change. At phase transition the speed of 
variation of these (fa) and (fa) and even their interrela- 
tion are changed. Some examples of this type gives toy 
model, considered in Appendix. 

Each of these transitions is accompanied by formation 
of bubbles of the old phase within the new one. Since 
cooling of Universe is very fast, these bubbles can be 
frozen in the new phase for a relatively long time (like 
in supercooled vapor). We can observe now some effects 
from different series of these bubbles despite the fact that 
these bubbles have disappeared by today. 

3) It is very attractive to assume that this very al- 
most degenerate sCPv state is realized now, ex- 
plaining observed CP violation. In this case at the first 
stage during the evolution of Universe domains of both 
types of the CP violation would form. Then, with mech- 
anisms like those discussed in [18] these domains decay 
to the phase with left violation of CP symmetry. In con- 
trast to the simple case considered in [18] domain wall 
can be high enough due to complex structure of poten- 
tial and the presence of a number of particles of matter 
within the walls. An important feature of the walls is 
that the matter within wall is much heavier than in vac- 
uum. Moreover, the profile of this wall in fa space can 
be complicated so that within the wall the photon can be 
massive (like in charged extremum) or not (if the saddle 
point between the minima obeys U(l) symmetry of elec- 
tromagnctism). In the former case the domain walls are 
opaque. In the second case these walls can be gray due 
to fluctuations. 

The comparison of speed of elimination of domains 
with the rate of cooling of Universe looks an interesting 
problem. 

In the considered case a small difference between the 
energies of two sCPv phases can have relation to the 
value of the cosmological constant. 
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APPENDIX A: APPENDIX. TOY MODEL 

To illustrate our general discussion we consider a sim- 
ple toy potential with weakly violated 4>\ <-> (j) 2 symme- 
try, where all the extrema can be calculated directly. In 
this model all coefficients are real and 

Ai = A 2 = A , A 3 = A(l - 5) , A 5 = Ax , 

m\ l =m 2 {l + A) , ml 2 = m 2 (l - A) , 

m 2 2 = m 2 xr/2 

at 5, A < 1 . 

In other words, our toy potential is 

v t = \ {xi+x 2 f - XSx lX2 + ^ (x 2 +xf) - 

TV? 

— ^(xi+x 2 )- 

m 2 A m 2 >cr 
7T- {xi-x 2 ) — {x 3 +x 3 ) . 



(Ala) 



Positivity constraint limits field of allowed parameters 
by inequality 



At S — A = this system has extra symmetry 
<Pi 4>2- This symmetry allows to obtain all extrema in 
analytical form but with specific degeneracies having no 
relation to reality. In the detailed analysis of solutions, 
considering A and 8 as perturbations, i.e. in the limit 
|A| ~ \S\ < 1, \A 2 /S\ < 1 we find that the values of the 
yi and extremum energies are changed by perturbation 
only weakly for neutral extrema. For charge extremum 
perturbations change values of yi strongly and diminish 
the range of possible realization of this extremum essen- 
tially while extremum energy shifts by perturbations only 
weakly. So that we present results of analysis for charged 
extremum at mentioned small but finite values of pertur- 
bations while for neutral extrema we consider solutions 
at A = 5 = only. 

It is useful in our analysis to use auxiliary quantities, 
which determine the scales of field and energy values at 
the extremum points, similar to SM 

F = m 2 /(2A), e = m 4 /(8A). 

• Charged extremum. 

First, we obtain values yi and energy from eqs. (20) 
and (22) in the first nontrivial approximation in A, 5: 



2/1,2 

£ch = 



v 

-e I 1 + 



r 

Y -, 

2' 

<5 2 +4A 2 ' 



2/3 



28 



(A2) 



At 5 = A = the additional tfii <-» cf> 2 symmetry results 
in degeneracy in values of y\ and y 2 , one can obtain only 
that yi + y 2 = Y . One can see that the values of fields 
at the extremum points depend strongly on the ratio of 
small quantities A/S while extremum energy depends on 
perturbations only weakly. 

The allowed values of yi i2 are limited by inequalities 
(6). These inequalities limit range of variables where this 
extremum can be realized by inequality 



1 > 



4A 2 



(A3) 



Note that this range determined for extra symmetry 
case A = 5 = is reduced strongly if A/S is not small. 
At 2A > 5 the charged extremum cannot be realized even 
if A and 5 are small. 

• sCPv extremum (d — A = 0). 

Using eqs. (33) we obtain j/i j2 ,cos£ and energy for 
sCPv extrema: 



Y 



Vi = 2/2 



2-x' 



cos£ 



£cpv 



+ 



(A4) 



2 > \h\ +S. 



(Alb) 



This extremum exists only when | cos£| < 1, i.e. 

|r(2-x)| <4. (A5) 
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Higgs masses. The mass of charged Higgs boson is 
given by eq. (31): 



M 2 



m 2 x 
2-x ' 



(A6) 



The neutral Higgs masses are eigenvalues of the mass 
matrix which we obtain using eq. (31): 



M = 



m 



A A - 



2-x 



A- . 
R 



A 



R 
R 



R 2(x+l)-2A (A7) 



with ,4 = l + ^cos 2 £, R = (xsin2£)/\/2. 

This form of matrix shows that one neutral Higgs 
is C-even, in terms of [4] that is hi = (rji — 172) /V~2 
(oei = 37r/4), its mass can be obtained by subtraction 
of the second line from the first one in the characteristic 
determinant: 



2 m 2 x 2 
M 2 Y = - = M 2 



2-x 



H± 



(A8a) 



Two other Higgs states have no definite CP parity they 
are mixed states of H = (771 -\-r\2)j\p2 and A with masses 



2(2 - x) 



(2 + x± 



±J (2 + xf - 8(2 - x)xsin 2 i 



(A8b) 



It is easy to check that the condition | cos£| < 1 (A5) 
at x > guarantee the positivity of all physical masses 
in model. 

• CPc extrema (8 = A = 0). 

To find yi's and energies for CPc extrema we rewrite 
eqs. (39) in the form (t = tan/3): 



(t 4 - 1) - 2t{t 2 - 1) = => 
(t 2 - l)[r(t 2 + 1) - 2t\ = , 



(A9) 



... 2\ 1 + rxt/2 
2 ' y ~ ' 1 'l + (l + x)t 2 



y="- = Y{l + t 2 ) ^ '- '"J 1 " = yi (l + 1 2 ) . (A10) 



Now the equation for extremum energy (38) takes form 

1 + rxt/2 



£cpc = -e{l+t 2 + xrt) 



1 + (1 + x)t 2 ' 



(All) 



The polynomial in the right-hand of (A9) is the prod- 
uct of two factors. The condition that one of this factors 
is describes different types of solution. 

Type I solutions are given by equation t 2 — 1 = 
and the extremum energy (All) is 

(J±): t = ±l^ yi =y 2 = Y 1 + r * K/2 , (A12) 

z + x 

(l+rtx/2) 2 

£l± = £ 1 + >,2/2 • (A13) 



Higgs masses. In this case the charged Higgs bo- 
son and the CP-odd Higgs boson masses are separated 
(see [4] M 2 H± - v 2 (v - A 5 /2), M 2 = v 2 (v - A 5 ), 
m\ 2 = 2v\Viv) 



, ,o xm 2 (rt — 1) 
M «±= 2 + x 



_ xm?[rt(2 - x) - 4] 
Ma 2(2 + x) 



(Al4a) 



while the mass matrix for the CP-even states and corre- 
sponding masses are 



M 



A B 



A{2 + x) \B A 

A = 4 + t(4 + x)rx , 

B = x 2 r + At(x+ 1) 

,,, xm 2 (rt-\) „ 

M 2 , = = M 2 

2 + x 

txr 



(A14b) 



M? = m 2 I 1 



At x > the conditions rt > 1, txr + 2 > and 
rt(2 — x) > 4 are necessary to have all mass squared to 
be positive, and for the extremum to be the minimum of 
the potential. The latter condition is opposite to (A5). 
Therefore sCPv extremum cannot exist if CPc extremum 
of type I is a minimum of potential. And vice versa if 
sCPv extremum exists, the CPc extremum of type I can- 
not be a minimum of potential. 

At x < the conditions rt < 1, txr + 2 > and 
rt(2 — x) < 4 are necessary to have all mass squared to 
be positive, and for the extremum to be the minimum of 
the potential. The latter condition coincides with (A5). 
Therefore if CPc extremum of type I is a minimum of 
potential, the sCPv extremum exists, but it is not a min- 
imum. 

Type II solutions break spontaneously (f>i <-> fa 
symmetry, they are given by second factor in (A9), i.e. by 
equation r(t 2 + 1) — 2t = 0, the solutions and extremum 
energy (All) are 



y h2 = I (1 ± VT^) , (A15) 



r 2 
£11 = £ C h (1 + C(A, S)) = 

^y) (1 + 0(A,J)) 



(A16) 



Second solution for the type (II) vacuum is obtained 
from that written here by the change of notations 
fa <-» fa- We do not consider it here. The real differ- 
ence appears in more realistic model violating fa <-> fa 
symmetry. 

Higgs masses. In this case the charged Higgs boson 
and the CP-odd Higss boson masses are separated too. 
Moreover, at S = A = we have M^ ± = due to 
accidental symmetry of model. So that we calculate this 
mass in the first order in small perturbations 5, A (for 
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other masses we neglect these perturbations): 
m 2 A 



M 2 H± 



M\ = - 



m 2 x 



For the CP even states we have 

kA m 2 (2 + x) (C r 
M = ■ 



C = l 

n 2 



4 \r C 
2-x 



2 + x 



(Al7a) 



(A17b) 



=* M%,h = X ( 2 + X ± 2)2 + 8 ^ 2 ) 

To have positive Mj| one must have x < while to 



have positive one must have x > 0. Therefore this 
cxtrcmum cannot be the minimum of the potential. 

• Case of negative m 2 . 

Up to this point we considered case to 2 > 0. However, 



the case 



< must be considered too. In this case 



Y < 0, and the conditions (6) allows only type (I) CP 
conserving extremum to exist. The constraints (6) for 
this case are: 



rtx<-2 (f = ±l). 



(A18) 



Otherwise the EW symmetry conserving solution 
(4>i) = is realized. 




2mf 2 /nT ! 



1B111P 

Neutral (I -) 




EW Symmery 








FIG. 2: Vacuum states in the plane xr = 2m? 2 /m 2 (vertical axis) — x — A 5 /A (horizontal axis). Left plot: m 2 > 0, right plot: 
m 2 < 0. 



• Comparison of extrema. Vacuum state. 

Taking into account the positivity constraint (Alb), we 
note that the condition r 2 < 1 is necessary for realization 
of the charged extremum and type II CPc extremum. It 
is sufficient for realization of the sCPv extremum as well. 

At different parameters of potential different ex- 
trema realize the vacuum state (with lowest energy) at 
A = S = 0. In particular, 



£cpv—£i = —£> { 



£ C h-£cPv = £ x — ' 



4(4 



2->c 



(A19) 



In the region where the charged extremum or CPc 
extremum II can exist, the sCPv extremum exists 
too. In this regions we have £ c h — Scpv > and 
£n — £cpv > 0. Therefore in the toy model the charged 



extremum and type II CPc extremum cannot be vacuum 
states. 

Now one must consider interrelation between the sCPv 
extremum and type I CPc extrema. 

According to (A19), at h < the difference 
£cpv — £i > 0, so the CPc extremum type I is realized in 
this region. On the other hand, at x > the difference 
£cpv — £i < 0, but the region where sCPv extremum 
exists is also constrained by eq. (A5). 

Fig. 2 shows regions with different vacuum states on 
the (x, nr) plane. The left plot shows different vacua for 
m 2 > 0. The right plot shows vacuum states for m 2 < 
0. In this case only type I CPc conserving extremum 
exists in the region constrained by (A18). In the region 
—2 < r x < 2, the EW symmetry conserving point yi = 
is the only minimum of the potential. 



